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Abstract 

Background: Understanding how human cells in tissue culture adapt to hypothermia may aid in developing new 
clinical procedures for improved ischemic and hypothermic protection. Human coronary artery endothelial cells 
grown to confluence at 37°C and then transferred to 25°C become resistant over time to oxidative stress and injury 
induced by 0°C storage and rewarming. This protection correlates with an increase in intracellular glutathione at 
25°C. To help understand the molecular basis of endothelial cold-adaptation, isolated proteins from cold-adapted 
(25°C/72 h) and pre-adapted cells were analyzed by quantitative proteomic methods and differentially expressed 
proteins were categorized using the DAVID Bioinformatics Resource. 

Results: Cells adapted to 25°C expressed changes in the abundance of 219 unique proteins representing a broad 
range of categories such as translation, glycolysis, biosynthetic (anabolic) processes, NAD, cytoskeletal organization, 
RNA processing, oxidoreductase activity, response-to-stress and cell redox homeostasis. The number of proteins 
that decreased significantly with cold-adaptation exceeded the number that increased by 2:1. Almost half of the 
decreases were associated with protein metabolic processes and a third were related to anabolic processes 
including protein, DNA and fatty acid synthesis. Changes consistent with the suppression of cytoskeletal dynamics 
provided further evidence that cold-adapted cells are in an energy conserving state. Among the specific changes 
were increases in the abundance and activity of redox proteins glutathione S-transferase, thioredoxin and 
thioredoxin reductase, which correlated with a decrease in oxidative stress, an increase in protein glutathionylation, 
and a recovery of reduced protein thiols during rewarming from 0°C Increases in S-adenosylhomocysteine 
hydrolase and nicotinamide phosphoribosyltransferase implicate a central role for the methionine-cysteine 
transulfuration pathway in increasing glutathione levels and the NAD salvage pathway in increasing the reducing 
capacity of cold-adapted cells. 

Conclusions: Endothelial adaptation to mild-moderate hypothermia down-regulates anabolic processes and 
increases the reducing capacity of cells to enhance their resistance to oxidation and injury associated with 0°C 
storage and rewarming. Inducing these characteristics in a clinical setting could potentially limit the damaging 
effects of energy insufficiency due to ischemia and prevent the disruption of integrated metabolism at low 
temperatures. 




Genomics 



Background 

Hypothermia is utilized in several clinical situations. 
Moderate-to-deep hypothermia (32°-15°C) is the primary 
method for delaying ischemic injury in patients during 
cardiovascular [1,2], neurovascular [3] or trauma surgery 
[4]. More severe hypothermia (< 10°C) is used to protect 
isolated tissues and organs for transplantation. Therapeu- 
tic mild hypothermia (34°-32°C) improves the neurologi- 
cal recovery of patients following a cardiac arrest [5,6], 
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stroke [7], or traumatic brain injury [8]. In each situation, 
the reduction in temperature decreases the metabolic 
demand for energy and slows the progression of injury. 
However, hypothermia also disrupts metabolic integra- 
tion and impairs important functional processes, particu- 
larly when it is severe. For example, cold-induced 
vascular injury is characterized by a loss of endothelial 
cell-cell contact [9,10], a release of inflammatory cyto- 
kines [11,12], an increase in the expression of adhesion 
molecules [11,12], impaired vasoactivity [13,14] and cell 
death [15]. Unfortunately, both the protective and harm- 
ful mechanisms induced by hypothermia vary with 
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temperature and are poorly understood and therefore its 
clinical utility remains limited. 

Many organisms, from prokaryotes to plants and ani- 
mals adapt to a range of temperatures and remain viable. 
Poikilotherms and mammalian hibernators adapt to cold 
and become resistant to injury from more prolonged and 
severe hypothermia and ischemia. Interestingly, the cells 
of non-hibernating mammals adapt to mild-moderate 
hypothermia in tissue culture [16-19] in what is likely a 
conserved response to cold and several mechanisms of 
cold-induced changes in gene expression have been identi- 
fied. These include a generalized inhibition of transcrip- 
tion and translation, an increase in transcription of some 
RNA-binding chaperones, alternative splicing of pre- 
mRNAs and preferential translation of mRNAs that have 
cold-inducible internal ribosome entry sites (IRESs) [19]. 
In our studies of cold-adaptation, human coronary artery 
endothelial cells (HCAECs) cultured at 25°C become pro- 
gressively more resistant over time to 0°C-injury and in 
particular to the oxidative stress induced by exposure to 
0°C and rewarming [16]. The molecular basis of the adap- 
tation remains largely unknown but the resulting protec- 
tion at 0°C is due, in part, to the sequestration of 
catalytically active iron [16]. The protection may also be 
associated with an increase in intracellular glutathione, an 
important antioxidant and signaling molecule of the cell, 
at 25°C. Glutathione (GSH) reacts directly with free radi- 
cals, participates in the reductive detoxification of hydro- 
gen peroxide and organic peroxides, serves as a co-factor 
in the enzymatic breakdown of xenobiotics and reacts 
with protein thiols to form mixed disulfides (P-SSG) 
under conditions of mild oxidative stress [20] . Protein glu- 
tathionylation is a reversible modification that provides a 
mechanism for protecting proteins from irreversible oxida- 
tive damage and for regulating protein function and 
thereby many diverse cellular processes such as cytoskele- 
tal organization, ionic homeostasis and the expression of 
genes involved in antioxidant defenses [21]. Understanding 
the process of cold-adaptation and the characteristics of 
the protected state in human cells may be beneficial to the 
development of new clinical strategies for improved 
ischemic and hypothermic protection. 

The aim of the present study was to characterize the 
cold-adapted state and the potential role of GSH in 
endothelial protection at 0°C. We performed quantitative 
proteomic analyses of human coronary artery endothelial 
cells collected before and after 72 h of culture at 25°C and 
classified the changes in protein abundance using the 
DAVID Bio informatics Resource [22,23]. The most signifi- 
cant biological processes associated with cold-adaptation 
were translation, glycolysis, RNA processing, actin fila- 
ment-based processes, and mRNA metabolic processes. 
Evidence of a regulated decrease in energy-consuming 
processes such as protein, DNA and fatty acid synthesis 



and suppression of cytoskeletal dynamics suggest that the 
cold-adapted state is hypometabolic. There was also evi- 
dence of a complementary increase in antioxidant protec- 
tion as there were significant changes in categories 
relevant to GSH such as oxidoreductase activity and cell 
redox homeostasis. Among the many changes in cold- 
adapted endothelial cells were increases in the abundance 
of thioredoxin, thioredoxin reductase and glutathione 
S-transferase, which correlated with increases in their 
respective activities, increases in protein glutathionylation, 
recovery of reduced protein thiols and diminished oxida- 
tive stress following rewarming from 0°C. 

Results 

Analysis of changes in protein expression at 25°C 

In order to characterize further the array of proteins 
altered during endothelial adaptation to cold, we used a 
label-free liquid chromatography/mass spectrometry (LC/ 
MS)-based method [24] to study protein extracts from 
HCAECs cultured at 37°C or cells subsequently cold- 
adapted to 25°C for 72 h as this length of exposure to 
25°C markedly increases GSH levels in cells and protects 
from 0°C cold exposure [16]. A total of 1089 proteins were 
identified (summary: Table 1; complete list of proteins: 
Additional file 1) and the number of proteins that 
decreased in abundance exceeded the number that 
increased by about 2:1 (202 decreased:102 increased). 
Priority 1 or 2 proteins were selected for further analysis 
and 181 of the 219 protein identifiers were successfully 
converted from the International Protein Index [25] to 
UniProt accession numbers. The assignment of protein 
priority is based on the peptide identification confidence 
from MS/MS data, which was described in detail pre- 
viously [24,26]. 

To aid in the interpretation of the experimental findings, 
the accession numbers were uploaded to the DAVID 
Bioinformatics Resources website [22,23] and proteins 
were categorized by Gene Ontology (GO) biological pro- 
cess, molecular function or cellular component terms [27] . 
Twenty-four representative categories of interest are 
shown in Table 2 (see Additional file 2 for the complete 
set of categories). A large number of proteins that were 
up- or down-regulated by hypothermia are involved in pri- 
mary metabolism (anabolic and catabolic processes), 
macromolecular metabolism, gene expression, biosynthesis 
(anabolic processes), protein metabolic processes and 
translation (56.7, 48.9, 34.8, 33.7, 31.5 and 22.5% of the 
total number of proteins that were analyzed, respectively). 
For these processes, the number of proteins that decreased 
exceeded the number that increased by ratios of 19-to-l 
(translation), 8.4-to-l (protein metabolism), 4.2-to-l 
(gene expression), 3.6-to-l (biosynthetic processes), 3-to-l 
(macromolecular metabolism), and 2-to-l (primary 
metabolic process). In addition, the categories of 
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Table 1 Changes in protein expression following cold-adaptation as determined by LC/MS-based label-free protein 



quantitative analysis 



Protein 
priority 


Peptide ID 
confidence 


Multiple 
sequences 


Number of 
proteins 


Number 
significant 
changes 


Maximum absolute 
fold change 


Number 
increasing at 
25°C 


Number 
decreasing at 
25°C 


1 


High 


Yes 


440 


169 (15.5%) 


1.70 


73 


96 


2 


High 


No 


364 


50 (4.6%) 


1.92 


12 


38 


3 


Moderate 


Yes 


7 


4 (0.04%) 


1.23 


0 




A 


Moderate 


No 


278 


81 (7.4%) 


2.37 


17 


6-1 


Overall 






1089 


304 (27.9%) 


2.37 


102 (9.4%) 


202 (18.5%) 



A total of 1089 proteins were identified by LC/MS and 27.9% of these had significantly increased or decreased in abundance following 25°C/72 h. Priorities 1-4 
were assigned based on confidence in peptide identification (High and Moderate) and whether or not there are multiple amino acid sequences (Yes or No) 
identified. Please see references [24] and [26] for more details about how the protein priority is determined. 



macromolecular complex assembly, ribonucleoprotein 
complex biogenesis, RNA processing, messenger RNA 
metabolic process, intracellular transport, protein folding, 
unfolded protein binding, microtubule cytoskeleton and 
response to stress also have more proteins that decreased 
in abundance than increased. Some processes show a net 
increase in proteins at 25°C. For example, for the catabolic 
process, the number of proteins that increased exceed the 



number that decreased by a ratio of 2.5-to-l. Twelve of 
the 21 catabolic proteins are also classified as glycolytic 
proteins, of which 1 1 increased and 1 decreased signifi- 
cantly at 25°C (Figure 1). GO categories such as cell prolif- 
eration, anti-apoptosis, cytoskeletal organization, actin 
filament-based processes, oxidoreductase activity and cell 
redox homeostasis also contain a greater number of pro- 
tein increases than decreases. 



Table 2 Selected GO a biological process classifications of proteins differentially expressed in HCAECs following 
cold-adaptation 



# of Proteins 



UP 


DN 


Ratio b 


% c 


Gene ontology categories 


P-Value d 


Fold 6 


34 


67 


-2.0 


56.7 


Primary metabolic process 


1 .4 X 1 0~ 3 


1.1 


22 


65 


-3.0 


48.9 


Macromolecule metabolic process 


1 .9 x 1 0~ 3 


1.3 


12 


50 


-4.2 


34.8 


Gene expression 


3.0 x 1 0~ 6 


1.7 


13 


47 


-3.6 


33.7 


Biosynthetic process 


1 .7 x 1 0~ 3 


1.4 


15 


6 


2.5 


11.8 


Catabolic process 


9.9 x 1 0~ 2 


1.4 


11 


1 


11.0 


6.7 


Glycolysis 


5.1 x 10~ 12 


21.6 


7 


59 


-8.4 


31.5 


Protein metabolic process 


3.9 x 1 0~ 5 


1.7 


2 


16 


-8.0 


10.1 


Macromolecular complex assembly 


2.2 X 1 0~ 3 


2.3 


0 


11 




6.2 


Ribonucleoprotein complex assembly 


5.3 X 10~ 5 


5.2 


9 


18 


-2.0 


15.2 


RNA processing 


1.1 X 10~ 9 


4.2 


8 


11 


-1.4 


10.7 


mRNA metabolic process 


3.6 X 1 0' 7 


4.3 


2 


38 


-19.0 


22.5 


Translation 


2.8 X 1 0~ 28 


10.2 


5 


11 


-2.2 


9.0 


Intracellular transport 


1.1 X 10~ 2 


2.1 


0 


7 




3.9 


Protein folding 


1 .8 X 1 0~ 2 


3.3 


0 


9 




5.1 


Unfolded protein binding' 


5.0 X 1 0~ 5 


6.9 


6 


5 


1.2 


6.2 


Cell proliferation 


3.3 X 10~ 2 


2.1 


6 


1 


6.0 


3.9 


Anti-apoptosis 


3.5 X 10~ 2 


2.9 


3 


3 


1.0 


3.4 


Structure-specific DNA binding 


2.5 X 1 0~ 2 


3.6 


1 1 


7 


1.6 


10.1 


Cytoskeletal organization 


1 .5 X 1 0~ 5 


3.5 


13 


4 


3.2 


9.6 


Actin filament-based process 


2.5 X 1 0~ 8 


6.0 


1 


11 


-11.0 


6.7 


Microtubule cytoskeleton 9 


3.3 X 1 0~ 2 


3.0 


9 


5 


1.8 


7.9 


Oxidoreductase activity' 


5.1 X 10~ 2 


1.8 


3 


2 


1.5 


2.8 


Cell redox homeostasis 


6.4 X 1 0~ 3 


6.7 


11 


16 


-1.5 


15.2 


Response-to-stress 


8.9 X 1 0~ 2 


1.4 



Proteins were classified using the Gene Functional Annotation Tool of DAVID bioinformatics resources, a. Gene Ontology; b. Ratios are negative where the 
number of proteins down-regulated at 25°C (# DN) > the number of proteins up-regulated at 25°C (# UP) and positive where # UP > # DN; c. % = the number of 
proteins in this category out of 178 proteins; d. Modified Fisher Exact P-value (probability of obtaining this many differentially expressed proteins in the category 
by chance); e. Fold-enrichment; f. GO Molecular Function term. g. GO Cellular Component term. 
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up-regulated by hypothermia 
down-regulated by hypothermia 



Gene 


UniProt 




Name 


Accession 


Protein Name 


ALDOA 


P04075 


aldolase A, fructose-bisphosphate 


ALDOC 


P09972 


aldolase C, fructose-bisphosphate 


EN01 


P06733 


enolase 1 , (alpha) 


EN02 


P09104 


enolase 2 (gamma, neuronal) 


EN03 


P13929 


enolase 3, (beta) 


GAPDH 


P04406 


glyceraldehyde-3-phosphate dehydrogenase 


LDHA 


P00338 


lactate dehydrogenase A 


LDHAL6A Q6ZMR3 


lactate dehydrogenase A-like 6A 


LDHAL6B Q9BYZ2 


lactate dehydrogenase A-like 6B 


PFKL 


P17858 


phosphofructokinase, liver 


PGK1 


P00558 


phosphoglycerate kinase 1 


PKM2 


P14618 


pyruvate kinase, muscle 



Figure 1 Glycolytic proteins that were significantly up- or down-regulated at 25°C. 



Two parameters used in DAVID, P-value and fold- 
enrichment, quantify the extent to which a given cate- 
gory is overrepresented within the dataset when com- 
pared to the precentage of similarly categorized genes 
from the entire genome. Categories with smaller P- 
values (< 0.05) and larger fold-enrichments (> 1.5) are 
considered to be of greater interest [23]. The most 
enriched groups from Table 2 are translation, glycolysis, 
RNA processing, actin filament-based processes and 
mRNA metabolic processes (P-values of 2.8 x 10' , 5.1 
x 10~ 12 , 1.1 x 10" 9 , 2.5 x 10* and 3.6 x 10" 7 , and fold- 
enrichment of 10.2, 21.6, 4.2, 6.0, and 4.3, respectively; 
see Table 2). Categories that may be relevant to GSH 
activity, such as oxidoreductase activity, response to 
stress and cell redox homeostasis appeared to be less 
enriched and of lower significance (fold-enrichment of 
1.8, 1.4, and 6.7 and P-values of 5.1 x 10~ 2 , 8.9 x 10~ 2 
and 6.4 x 10' 3 , respectively) and therefore potentially of 
reduced importance. These categories show mixed 
changes in protein abundance, with increase-to-decrease 
ratios of 1.8-to-l, l-to-1.5, and 1.5-to-l. 

Oxidative stress and protein oxidation/reduction 

Figure 2 lists proteins that have oxidoreductase activity 
or regulate cell redox homeostasis and were significantly 



up- or down-regulated at 25°C. Thioredoxin (TRX1), 
thioredoxin reductase 1 (TXNRD1), peroxiredoxin-1 
(PRDX1) and glutathione S-transferase (GST, omega 1) 
increased in abundance and protein disulfide isomerases 
1 and 4 (PDIA1, PDIA4) decreased significantly. TRX1 
is an important cytoplasmic protein that reduces oxi- 
dized proteins such as peroxiredoxins and ribonucleo- 
tide reductase [28]. Thioredoxins (Trx) are reduced by 
thioredoxin reductases (TrxR) with reducing equivalents 
from NADPH in a manner analogous to the reduction 
of oxidized to reduced glutathione by GSSG reductase 
and NADPH. Peroxiredoxins play a protective role in 
cells by eliminating peroxides [28] and GSTs are impor- 
tant for detoxifying the products of oxidative stress, 
such as oxidized proteins, DNA or lipids, through their 
conjugation with GSH [20]. PDIA's 1 and 4 are found in 
the endoplasmic reticulum (ER) and are involved in the 
oxidative folding of nascent proteins by inducing the 
formation of intramolecular disulfide bonds [29]. The 
proteomic changes indicate that hypothermic adaptation 
potentially increases the cell's capacity to reduce oxi- 
dized proteins. 

To determine if cold-adapted cells have an increased 
capacity to reduce oxidized proteins, confluent endothe- 
lial cells were cultured at 25°C for 72 h or maintained at 
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up-regulated by hypothermia 
down-regulated by hypothermia 



Gene 


UniProt 




Name 


Accession 


Protein Name 


SAHH 


P23526 


s-adenosylhomocysteine hydrolase 


IMPDH2 


P12268 


IMP (inosine monophosphate) dehydrogenase 2 


ALDH1A1 


P00352 


aldehyde dehydrogenase 1 family, member A1 


FAS 


P49327 


fatty acid synthase 


GST01 


P78417 


glutathione s-transferase omega 1 


GAPDH 


P04406 


glyceraldehyde-3-phosphate dehydrogenase 


LDHA 


P00338 


lactate dehydrogenase A 


LDHAL6A 


Q6ZMR3 


lactate dehydrogenase A-like 6A 


LDHAL6B 


Q9BYZ2 


lactate dehydrogenase A-like 6B 


PRDX1 


Q06830 


peroxiredoxin 1 


PLOD2 


000469 


procollagen-lysine, 2-oxoglutarate 5-dioxygenase 2 


PDIA1 


P07237 


protein disulfide isomerase 


RRM1 


P23921 


ribonucleotide reductase ml polypeptide 


TXNRD1 


Q16881 


thioredoxin reductase 1 


YWHAZ 


P63104 


14-3-3 protein zeta 


PDIA4 


P13667 


protein disulfide isomerase family 1 , member 4 


TRX1 


P10599 


thioredoxin 



Figure 2 Proteins that participate in redox reactions and were significantly up- or down- regulated at 25°C 



37°C (controls) and protein thiols (PSH), protein glu- 
tathionylation (PSSG) and free glutathione (oxidized, 
GSSG and total, tGSH) was measured immediately 
before or after exposure of cells to 0°C and subsequent 
rewarming. HCAECs are strongly susceptible to oxida- 
tive injury induced by 0°C exposure and rewarming [16]. 
PSH and PSSG in cells was unchanged after 72 h at 25° 
C (Figures 3A, B) as the treatment did not induce 
changes in GSSG/tGSH (Figure 3C). Exposure of control 
cells to 0°C followed by rewarming to 37°C for 3 h, 
however, decreased PSH (*P < 0.01; Figure 3A) and 
increased GSSG/tGSH (*P < 0.001 Figure 3C), which is 



indicative of oxidative stress. Similarly, cells treated with 
the pro-oxidant diamide at 37°C (1 mM/30 min) showed 
a significant decrease in PSH and an increase in GSSG/ 
tGSH (data not shown). However, PSH in cold-adapted 
cells was significantly greater than in non-adapted cells 
following rewarming from 0°C (**P < 0.001 vs. 37°C+0°C 
+RW; Figure 3A). This may be indicative of enzymatic 
reduction of oxidized protein thiols during rewarming 
as PSH increased significantly during this phase (+P < 
0.05 vs. 25°C+0°C; Figure 3A). PSSG trended lower with 
0°C storage and 0°C plus rewarming but not in cold- 
adapted cells, which trended oppositely (*P < 0.05 vs. 



Zieger ef al. BMC Genomics 201 1, 12:630 
http://www.biomedcentral.eom/1 471 -2 1 64/1 2/630 



Page 6 of 22 




control 0°C/48 h 0°C+RW 




control 0°C/48 h + RW 



Figure 3 Effect of cold-adaptation (25°C/72 h) and cold storage 

(0°C/48 h) on protein and non-protein thiols (A) Cold-adaptation 
prevents the significant loss of protein thiols (PSH) that occurs 
during cold storage/rewarming (*P < 0.01 vs. 37°C control; ** P < 
0.001 vs. 37°C+0°C+RW; n = 5 experiments), possibly by the 
reduction of oxidized proteins during the rewarming phase (+P < 
0.05 vs. 25°C+0°C; n = 5 experiments). (B) Cold-adapted cells have a 
significantly higher level of glutathiolated protein (PSSG) than non- 
adapted cells (37°C) following cold storage/rewarming (**P < 0.05 
vs. 37°C+0°C+RW; n = 5 experiments). (C). Cold-adaptation 
attenuates the increase in GSSG/tGSH, indicative of oxidative stress, 
that occurs during cold storage/rewarming (*P < 0.001 vs. 37°C 
control; **P < 0.001 vs. 37°C+0°C+RW; n = 3 experiments). Open 
bars represent cell cold-adapted at 25°C for 72 h; closed bars 
represent cells maintained at 37°C for 72 h. 



37°C+0°C+RW; Figure 3B). Thus, cold-adapted cells 
have an increased capacity to recover PSH following 
cold-storage and rewarming when compared to non- 
adapted cells as well as an increased ability to glutathio- 
nylate proteins, which correlates with the observed 
increase in GST. 

To determine if the Trx-TrxR and GSH-GST systems 
are functionally more active in cold-adapted endothelial 
cells, Trx, TrxR and GST activities were measured in cells 
adapted to 25°C for 72 h and in 37°C control cells before 
and following 0°C storage plus rewarming. Figure 4A 
shows there is a significant increase in Trx activity with 
cold-adaptation (*P < 0.05 vs. 37°C control) and a signifi- 
cantly higher level of Trx and TrxR activity following 
rewarming from 0°C when compared to cold-stored con- 
trol cells (**P < 0.01 vs. 37°C + 0°C/48 h + RW, Figure 4A 
and 5! P < 0.05 vs. 37°C + 0°C/48 h + RW, Figure 4B). This 
correlates with the increased recovery of PSH in cold- 
adapted cells following rewarming from 0°C (Figure 3A). 
GSH-dependent GST activity is also increased in cold- 
adapted cells following 0°C storage and rewarming when 
compared to cold-stored non-adapted cells ("P < 0.05 vs. 
37°C + 0°C/48 h + RW; Figure 4C), which correlates with 
the increase in PSSG (Figure 3B) and tGSH observed fol- 
lowing rewarming from 0°C [16]. In summary, the findings 
indicate that the increase in Trx-TrxR and GSH-GST 
activities may be an important component of the antioxi- 
dant protection of cold-adapted endothelial cells and their 
proteins during 0°C storage and rewarming. 

NAD metabolic proteins 

The molecule nicotinamide adenine dinucleotide or 
NAD plays a central role in linking the energy produ- 
cing catabolic reactions of the cell with the maintenance 
of redox homeostasis [30]. Significant changes in NAD 
metabolic proteins were identified in this study. Cells 
adapted to 25°C had an increase in the level of nicotina- 
mide phosphoribosyltransferase (Nampt), the rate-limit- 
ing enzyme in NAD biosynthesis from nicotinamide [31] 
(see Additional file 1). There was also an increase in gly- 
colytic enzymes that use NAD as an electron acceptor 
or cofactor. Glyceraldehyde-3-phosphate dehydrogenase 
(GAPDH), which reduces NAD + to NADH, increased at 
25°C, as did lactate dehydrogenase (LDHA, LDHAL6A 
and LDHAL6B), which reduces pyruvate to lactate and 
regenerates NAD + (Figure 2). Interestingly, cold-adapta- 
tion increased the abundance of S-adenosylhomocys- 
teine hydrolase (SAHH), an NAD + dependent enzyme 
and component of the methionine-cysteine transulfura- 
tion pathway that generates cysteine used to synthesize 
GSH [32]. Inosine monophosphate dehydrogenase 
(IMPDH2), the NAD + -dependent and rate-limiting 
enzyme of de novo guanine nucleotide synthesis [33], 
decreased at 25°C. 
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Figure 4 Effect of cold-adaptation (25°C/72 h) and cold storage 
(0°C/48 h) on thioredoxin (Trx), thioredoxin reductase (TrxR) 
and GST (glutathione S transferase) activity (A) Cold-adaptation 
increases Trx activity before or after 0°C storage (*P < 0.05 vs. 37°C 
control; **P < 0.001 vs. 37°C + 0°C/48 h + RW; n = 5 experiments). 
(B) TrxR activity is greater in cold-adapted cells following 0°C 
storage and rewarming (*P < 0.05 vs. 37°C + 0°C/48 h + RW; n = 4 
experiments). (C) GST activity in cold-adapted cells following 0°C 
storage and rewarming exceeds that of control cells (*P < 0.05 vs. 
37°C + 0°C/48 h + RW; n = 4 experiments). Open bars represent cell 
cold-adapted at 25°C for 72 h; closed bars represent cells 
maintained at 37°C for 72 h. 



The phosphorylated and reduced form of NAD, 
NADPH, provides the reducing equivalents for GSH and 
thioredoxin reduction that bolsters the cell's resistance 
to oxidative stress. Several NADP(H)-dependent 
enzymes changed significantly in abundance during 
cold-adaptation. Aldehyde dehydrogenase (ALDH1A1), 
an NADP + -dependent enzyme, increased at 25°C (Figure 
2). Its function in the cell is to oxidize toxic aldehydes 
produced during lipid peroxidation and thereby prevent 
their adduct formation with proteins or DNA [34]. Con- 
versely, several proteins decreased in abundance, includ- 
ing ribonucleotide reductase (RRM1), a thioredoxin- 
dependent enzyme [35] that catalyzes the formation of 
deoxyribonucleotide from ribonucleotide [36]; and fatty 
acid synthase (FAS), which catalyzes the formation of 
long-chain fatty acids from acetyl CoA, malonyl CoA 
and NADPH. 

Proteins that regulate gene expression 

Cold-adaptation induced significant changes in the abun- 
dance of proteins that organize chromatin, bind to RNA 
and DNA, and regulate transcription and translation 
(Figure 5). There were changes in the abundance of histone 
proteins (H2AV, H2BFC, H2BFQ), non-histone chromatin- 
binding proteins (HMGB1, HMG1L2, BANF1) and the his- 
tone chaperone NASP at 25°C Two proteins that increased 
in abundance at 25°C induce transcriptional repression 
(ENOl and HNRNPDL) and 2 cold shock domain proteins 
(YB-1, CSDA), which repress growth factor and stress 
response genes [37,38], decreased at 25°C. Another tran- 
scription factor that decreased at 25°C, |3-catenin, binds to 
and regulates ROS-sensitive transcription factors such as 
HIF-1 and FOXO [39]. There was also a decrease in 11 
ribosomal proteins and proteins associated with transcrip- 
tion initiation (BTF3), DNA replication/repair (PCNA), 
pre-mRNA processing (NONO, SRFS1, HNRNPH1, 
HNRNPF, CPSF6), mRNA elongation (HNRNPU), non- 
sense-mediated mRNA decay (EIF4A3), stress granule 
assembly (G3BP1, G3BP2) [40], translation regulation 
(RPS14, PABPC4, PABPC1) and targeting of nascent pro- 
teins to the endoplasmic reticulum (NACA). Conversely, 
there was a significant increase in RNA-binding proteins 
(PTBP1, HNRNPL and ADARB1), which process pre- 
mRNA, and HNRNPK, HNRNPA1 and PCBP2, which are 
multi-functional proteins involved in such processes as 
chromatin remodeling, transcriptional repression/activa- 
tion, mRNA export, stabilization and turnover, and transla- 
tional silencing/activation [41]. RBM3, a cold-inducible 
RNA-binding protein that is believed to enhance global 
protein synthesis [42], increased significantly at 25°C. 

Response to stress proteins 

Figure 6 lists 27 response-to-stress proteins that were sig- 
nificantly up- or down-regulated with cold-adaptation. 
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down -regulated by hypothermia 



Gene 


UniProt 




Name 


Accession 


Protein Name 


YB-1 


P67809 


y box binding protein 1 


CSDA 


P16989 


cold shock domain protein a 


EN01 


P06733 


enolase 1, (alpha) 


HNRNPDL 


014979 


heterogeneous nuclear ribonucleoprotein d-like 


HMG1L2 


P51858 


high-mobility group protein 1 -like 2 


RPS14 


P62263 


ribosomal protein s14 


CTNNB1 


P35222 


catenin (cadherin-associated protein ), beta 1, 88kda 


HMGB1 


P09429 


high-mobility group box 1 


NONO 


Q15233 


non-pou domain containing, octamer-binding 


FLNA 


P21333 


filamin a, alpha (actin binding protein 280) 


CALR 


P27797 


calreticulin 


BTF3 


P20290 


basic transcription factor 3 


YWHAQ 


P27348 


14-3-3 protein theta 


RHOQ 


P17081 


ras homolog gene family, member q 


NACA 


Q13765 


nascent-polypeptide-associated complex alpha polypeptide 


H2AV 


Q71UI9 


h2a histone family, member v 


H2BFC 


Q99880 


histone 1 , h2bl 


NASP 


P49321 


nuclear autoantigenic sperm protein (histone-binding) 


EIF4A3 


P38919 


eukaryotic translation initiation factor 4A, isoform 3 


G3BP1 


Q 13283 


ras-gtpase-activating protein sh3-domain binding protein 


H2BFQ 


Q16778 


histone 1 , h2bb 


PCNA 


P 12004 


proliferating cell nuclear antigen 


HNRNPU 


Q00839 


heterogeneous nuclear ribonucleoprotein u (scaffold attachment factor a) 


BANF1 


075531 


barrier to autointegration factor 1 


HNRNPA1 


P09651 


heterogeneous nuclear ribonucleoprotein a1 


HNRNPK 


P61978 


heterogeneous nuclear ribonucleoprotein k 


PCBP2 


Q 15366 


poly(rc) binding protein 2 


HNRNPH1 


P31943 


heterogeneous nuclear ribonucleoprotein hi 


PABPC4 


Q13310 


poly(a) binding protein, cytoplasmic 4 (inducible form) 


RBM3 


P98179 


rna binding motif (rnpl, rrm) protein 3 


PABPC1 


P11940 


poly(a) binding protein, cytoplasmic 2 


SFRS1 


Q07955 


splicing factor, arginine/serine-rich 1 (splicing factor 2, alternate splicing factor) 


HNRNPF 


P52597 


heterogeneous nuclear ribonucleoprotein f 


CPSF6 


Q 16630 


cleavage and polyadenylation specific factor 6, 68kd 


HNRNPL 


P14866 


heterogeneous nuclear ribonucleoprotein I 


PTBP1 


P26599 


polypyrimidine tract binding protein 1 


ADARB1 


P78563 


adenosine deaminase, rna-specific, b1 


RPL7 


P18124 


ribosomal protein L7 


RPS27 


P42677 


ribosomal protein S27 


RPS3 


P23396 


ribosomal protein S3 


RPL10A 


P62906 


ribosomal protein L10a 


RPS17 


P08708 


ribosomal protein S17 


RPS19 


P39019 


ribosomal protein S19 


RPS24 


P62847 


ribosomal protein S24 


RPS28 


P62857 


ribosomal protein S28 


RPS6 


P62753 


ribosomal protein S6 


RPS7 


P62081 


ribosomal protein S7 



Figure 5 Proteins that may regulate gene expression and were significantly up- or down-regulated at 25°C. 



Eight are molecular chaperones or chaperone-associated 
proteins that assist with protein folding or macromolecu- 
lar assembly and all were down-regulated (HYOU1, 
HSPA7, HSPA1A, HSP90B, HSPA8, STIP1, DNAJA1, 



and HSP90AA2). Other proteins that play a role in pro- 
tein folding decreased in abundance, including PDIA1 
and PDIA4, chaperonins CCT4 and CCT5, the co-cha- 
perone CDC37 and the calcium-binding chaperone 
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B up-regulated by hypothermia 
down-regulated by hypothermia 



Gene 


UniProt 




Name 


Accession 


Protein Name 


TXNRD1 


Q16881 


thioredoxin reductase 1 


PDIA4 


P13667 


protein disulfide isomerase family 1, member 4 


PDIA1 


P07237 


protein disulfide isomerase 


PLOD2 


000469 


procollagen-lysine, 2-oxoglutarate 5-dioxygenase 2 


PRDX1 


Q06830 


peroxiredoxin 1 


YWHAZ 


P63104 


14-3-3 zeta 


FN1 


P02751 


fibronectin 1 


NONO 


Q 15233 


non-pou domain containing, octamer-binding 


HYOU1 


Q9Y4L1 


hypoxia up-regulated 1 


PABPC4 


Q13310 


poly(a) binding protein, cytoplasmic 4 (inducible form) 


HSPA7 


P48741 


heat shock 70kda protein 1 -like 


CTGF 


P29279 


connective tissue growth factor 


HSPA1A 


P08107 


heat shock 70kda protein 1a 


PAI1 


P05121 


serpine peptidase inhibitor, clade e, member 1 


CSDA 


P 16989 


cold shock domain protein a 


HMGB1 


P09429 


high-mobility group box 1 


PCNA 


P12004 


proliferating cell nuclear antigen 


HSP90B 


P08238 


heat shock protein 90kda alpha (cytosolic), class b member 1 


HSPA8 


P11142 


heat shock 70kda protein 8 


STIP1 


P31948 


stress-induced-phosphoprotein 1 (hsp70/hsp90-organizing protein) 


TSP1 


P07996 


thrombospondin 1 


DNAJA1 


P31689 


dnaj (hsp40) homolog, subfamily a, member 1 


HSP90AA2 


Q 14568 


heat shock protein 90kda alpha (cytosolic), class a member 2 


ANXA1 


P04083 


annexin a1 


CCT4 


P50991 


chaperonin containing tcpl, subunit 4 (delta) 


CALR 


P27797 


calreticulin 


CCT5 


P48643 


chaperonin containing tcpl, subunit 5 (epsilon) 


CDC37 


Q 16543 


hsp90 co-chaperone cdc37 


PAPSS2 


095340 


3'-phosphoadenosine 5'-phosphosulfate synthase 2 


RBM3 


P98179 


RNA binding motif protein 3 


ALD03 


P09972 


aldolase C, fructose-bisphosphate 


EN03 


P 13929 


enolase 3 (beta, muscle) 


H2BFQ 


Q 16778 


histone cluster 2, H2be 


RPS3 


P23396 


ribosomal protein S3 



Figure 6 

Activity is 



Response-to-stress 

in Figure 2). 



proteins that were significantly up- or down-regulated at 25°C (Note: the complete list for Oxidoreductase 



CALR. Five of the proteins listed in Figure 6 are ER-resi- 
dent proteins (PDIA1, PDIA4, PLOD2, HYOU1, CALR) 
and all decreased at 25°C. Other stress response proteins 
that were down-regulated are associated with DNA 



replication/repair (PCNA), transcriptional repression 
(CSDA), pre-mRNA processing (NONO) and regulation 
of translation (PABPC4). There was, however, up-regula- 
tion of the non-histone chromatin-binding protein 
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HMGB1, the anti-inflammatory mediator annexin Al 
(ANXA1), the adapter protein 14-3-3 zeta and 3 proteins 
associated with cell adhesion (FN1, TSP1, CTGF). 

Cytoskeletal proteins 

The proteomic findings indicate that there were significant 
changes in the cytoskeletal organization of cold-adapted 
cells, including changes in the abundance of 36 proteins 
associated with microtubule, actin filament and intermedi- 
ate filament networks (Figure 7). Eleven out of 12 identi- 
fied proteins of the microtubule cytoskeleton decreased 
significantly, including a-tubulins al, a2, a3, a3E, and 
a4, as well as proteins involved in microtubule assembly 
(RANBP1), transport (CKAP5) and stability (CKAP2). The 
decrease in a-tubulin abundance was matched by a 
decrease in chaperonins CCT4 and CCT5 (Figure 6), 
which assist in actin and tubulin folding. 

Thirteen out of 17 proteins associated with actin fila- 
ment (F-actin) -based processes were up-regulated in cold- 
adapted cells (Figure 7). Although there were no signifi- 
cant changes in the abundance of cytosolic actin at 25°C 
(see Additional file 1), there were significant decreases in 
actin nucleation and cross-linking proteins (ARPC3 [43] 
and FLNA [44]), increases in actin filament severing pro- 
teins (CFL1 and DSTN [45]) and in proteins that inhibit 
polymerization by binding to monomeric G-actin (CAP1 
[46], TMSB10 [47] and TMSL3). There was also a 
decrease in septin 2, a GTP-binding and actin-organizing 
protein [48]. Cold- adaptation increased several other pro- 
teins that associate with actin, including the membrane- 
associated small GTPase RHOQ; Rho GDP dissociation 
inhibitor p (ARHGDIB), a chaperone protein that inacti- 
vates Rho GTPases [49]; CDC42, a Rho GTPase that regu- 
lates actin cytoskeletal remodeling; the actin-bundling 
protein fascin (FSCN1); and myosin light chain subunits 
(MYL6, MYL6B, MYL12A and MYL12B), which regulate 
actinomyosin contractility. 

Other changes included a cold-induced increase in the 
type III intermediate filament (IF) proteins vimentin 
(VIM) and peripherin (PRPH) that form a third class of 
cytoskeletal network (Figure 7). There were changes in the 
abundance of proteins (CASL, TLN1, PLEC1, and 
FERMT3) that link the cytoskeleton to integrin, which is a 
membrane receptor that binds the cell to proteins of the 
extracellular matrix, and there was a significant decrease 
in P-catenin (CTNNB1), a cytoplasmic component of the 
adherens cell-cell junction. 

Discussion 

Oxidative stress and protein oxidation/reduction 

Differential proteomic analyses of endothelial cells col- 
lected before or after cold-adaptation at 25°C demon- 
strated significant changes in the abundance of proteins 
that regulate redox homeostasis or have oxidoreductase 



activity. The increase in thioredoxin and thioredoxin 
reductase correlated with the increase in their respective 
activities and with the recovery of reduced protein thiols 
following rewarming from 0°C, conditions that generate 
oxidative stress and endothelial cell damage [16,50,51]. 
Protein disulfide isomerases catalyze thiol oxidation; 
therefore, the decrease in PDIA1 and PDIA4 abundance 
in cold-adapted cells may have contributed to the recov- 
ery of PSH with rewarming. Additional adaptations are 
potentially beneficial to the maintenance of redox home- 
ostasis. For example, preliminary studies showed that 
there is an increase in tyrosine hydroxylase, the rate-lim- 
iting enzyme in catecholamine synthesis, at 25°C [52]. 
Catecholamines such as dopamine attenuate hypother- 
mia-induced reactive oxygen species and the loss of -SH 
reducing equivalents at 4°C in human umbilical vein 
endothelial cells [53]. 

Endothelial cold-adaptation at 25°C is accompanied by a 
significant increase in intracellular GSH [16]. The two 
enzymes required for GSH synthesis, gamma glutamylcys- 
teine synthetase and glutathione synthetase, were not 
identified in our proteomic study. However, there was a 
significant increase in S-adenosylhomocysteine hydrolase 
(SAHH; Figure 2), a component of the methionine- 
cysteine transsulfuration pathway that generates cysteine 
for GSH synthesis. This pathway is up-regulated in 
response to oxidative stress [54] and accounts for up to 
50% of GSH production in some tissues [32] . We did not 
observe an increase in GSSG/tGSH, which is indicative of 
oxidative stress, with cold-adaptation, but this does not 
rule out the possibility of localized changes in GSSG/ 
tGSH that could have had a stimulatory effect on tGSH 
synthesis at 25°C. The increase in protein glutathionylation 
following rewarming from 0°C matched the higher levels 
of tGSH in cold-adapted cells [16] and the increase in 
abundance and activity of glutathione S-transferase. GST 
conjugates GSH to proteins under conditions of oxidative 
stress [55]. Protein glutathionylation is a reversible process 
that precludes irreversible disulfide cross-linking [56]; 
therefore, glutathionylation is potentially a protective 
mechanism in cold-adapted cells. 

NAD metabolic proteins 

The increase in GSH/GST and Trx/TrxR activities 
requires the availability of reducing equivalents from 
NADPH. Cold-adapted cells had a significant increase in 
Nampt, the rate-limiting enzyme in NAD synthesis from 
nicotinamide [30]. NAD kinase, which phosphorylates 
NAD + (H) to NADP + (H), was not identified in our study. 
The enzymes of the pentose phosphate pathway (PPP) 
that reduce NADP + to NADPH, glucose-6-phosphate 
dehydrogenase and 6-phosphogluconate dehydrogenase, 
were unchanged in abundance. However, the decrease in 
phosphofructokinase (PFK), the rate-limiting enzyme in 
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j up-regulated by hypothermia 
down-regulated by hypothermia 



Gene 


UniProt 




Name 


Accession 


Protein Name 


TUBA2 


Q 13748 


tubulin, alpha 2 


RANBP1 


P43487 


ran binding protein 1 


TUBA3E 


Q6PEY2 


similar to alpha tubulin 


CKAP5 


Q14008 


cytoskeleton associated protein 5 


TUBA3 


Q71U36 


tubulin, alpha 3 


TUBA4 


Q9H853 


tubulin, alpha 4 


TUBA1 


P68366 


tubulin, alpha 1 (testis specific) 


SEPT2 


Q15019 


septin 2 


CLIC5 


Q9NZA1 


chloride intracellular channel 5 


CTNNB1 


P35222 


catenin (cadherin-associated protein), beta 1, 88kda 


CASL 


Q14511 


neural precursor cell expressed, developmentally down-regulated 9 


DYNLL1 


P63167 


dynein, light chain, Ic8-type 1 


CAP1 


Q01518 


cap, adenylate cyclase-associated protein 1 (yeast) 


FSCN1 


Q16658 


fascin homolog 1, actin-bundling protein (stronglylocentrotus purpuratus) 


MYL6 


P60660 


myosin, light polypeptide 6, alkali, smooth muscle and non-muscle 


ARHGDIB 


P52566 


rho gdp dissociation inhibitor (gdi) beta 


CFL1 


P23528 


cofilin 1 (non-muscle) 


MYL6B 


P14649 


myosin light chain 1 slow a 


DSTN 


P60981 


destrin (actin depolymerizing factor) 


ARPC3 


015145 


actin related protein 2/3 complex, subunit 3, 21 kda 


FLNA 


P21333 


filamin a, alpha (actin binding protein 280) 


CALR 


P27797 


calreticulin 


CDC42 


P60953 


cell division cycle 42 (gtp binding protein, 25kda) 


TMSB10 


P63313 


thymosin, beta 10 


VIM 


P08670 


vimentin 


PRPH 


P41219 


peripherin 


PLEC1 


Q15149 


plectin 1, intermediate filament binding protein 500kda 


TLN1 


Q9Y490 


talin 1 


FERMT3 


Q86UX7 


fermitin family homolog 3 


MYL12A 


P19105 


myosin regulatory light chain mrcl3 


ANXA1 


P04083 


annexin a1 


ALDOA 


P04075 


aldolase A, fructose-bisphosphate 


ALDOC 


P09972 


aldolase C, fructose-bisphosphate 


HNRNPH1 


P31943 


heterogeneous nuclear ribonucleoprotein H 


MYL12B 


014950 


myosin regulatory light chain 12B 


RHOQ 


P17081 


ras homolog gene family, member q 


TMSL3 


A8MW06 


thymosin-like 3 


YWHAZ 


P63104 


14-3-3 zeta 



Figure 7 Cytoskeletal proteins that were significantly up- or down-reguiated at 25°C. 



the glycolytic pathway, potentially redirects glucose cata- 
bolism from glycolysis to the alternative PPP in a manner 
similarly induced by down-regulating the activity of other 
downstream glycolytic enzymes such as GAPDH or 



triosephosphate isomerase (TPI) [57]. The proteomic data 
indicate that there may be two additional mechanisms in 
cold-adapted cells to bolster NADPH availability. Alde- 
hyde dehydrogenase (ALDH1A1- Figure 2), which oxidizes 
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the lipid peroxidation products 4-hydroxynonenal and 
malondialdehyde using NADP + as a cofactor [34], 
increased at 25°C. Lipid peroxidation is an important 
mechanism of injury during 0°C storage and rewarming 
[50,51]. An increase in ALDH1A1 would therefore poten- 
tially protect cell membranes from oxidative damage and 
simultaneously regenerate NADPH for GSSG and Trx 
reduction. A second approach to ensuring NADPH avail- 
ability relates to the apparent decrease in NADPH-depen- 
dent biosyntheses at 25°C. Ribonucleotide reductase 
(RRM1), a Trx-dependent enzyme [35] that catalyzes 
deoxyribonucleotide formation and controls the rate of 
DNA synthesis at 37°C [36], decreased significantly at 25° 
C (Figure 2). Inosine monophosphate dehydrogenase 
(IMPDH2), the rate limiting enzyme in de novo guanine 
nucleotide synthesis, also decreased at 25°C. Inhibition of 
IMPDH2 decreases guanine nucleotide pools and inter- 
rupts both DNA and RNA synthesis [58]. Down-regulating 
these two enzymes may therefore suppress DNA synthesis 
independently of the Q 10 effects [59] and increase the 
availability of NADPH for GSSG and Trx reduction. A sig- 
nificant decrease in the NADPH-dependent fatty acid 
synthase (Figure 2), and presumably fatty acid synthesis, at 
25°C may have a similar NADPH-sparing effect and 
thereby contribute to the protective phenotype of cold- 
adapted cells. 

Glycolysis 

P-values and fold-enrichment values derived from the 
bioinformatic analysis (Table 2) indicate that one of the 
most important adaptations at 25°C may be changes to 
glycolysis. There are 10 steps in the conversion of glu- 
cose to pyruvate [60] and the enzymes that catalyze 6 of 
the final 8 steps changed significantly in abundance with 
cold-adaptation (Figure 8). Contrary to the aforemen- 
tioned decrease in PFK, however, there was an increase 
in the abundance of enzymes that catalyze 5 of the 
remaining 7 steps in glycolysis downstream of PFK, 
including aldolase, glyceraldehyde phosphate dehydro- 
genase, phosphoglycerate kinase, enolase, and pyruvate 
kinase (ALDOA/ALDOC, GAPDH, PGK1, ENOl/ 
EN02/EN03 and PKM2; Figure 2). This inconsistency 
is reconcilable if there is an increase in intermediate 
metabolites entering the glycolytic pathway. For exam- 
ple, glycerol is produced by the breakdown of triglycer- 
ides to fatty acids and is converted into glycerol 3- 
phosphate and then into dihydroxyacetone phosphate 
(DHAP), a glycolytic intermediate. The conversion of 
glycerol to DHAP and further oxidation to pyruvate 
would yield one mole of ATP per mole of glycerol con- 
sumed compared to two moles of ATP per mole of glu- 
cose. Glycerol is a substrate for gluconeogenesis under 
some conditions [61] but it has also been identified as 
an important energy substrate in the ex vivo rodent 



heart [62,63], cultured cardiomyocytes [64] and the 
rodent brain [65]. A role for glycerol in endothelial 
energy metabolism during adverse conditions has not 
been reported. 

Cold-adapted cells had a significant increase in lactate 
dehydrogenase (LDHA, LDHAL6A and LDHAL6B), 
which reduces pyruvate to lactate. This step diverts pyr- 
uvate from further oxidation by mitochondrial respira- 
tion and regenerates NAD + to facilitate more glycolysis. 
The reduction of pyruvate despite the availability of oxy- 
gen is known as the Warburg effect or aerobic glycolysis 
and has been shown to decrease mitochondrial respira- 
tion, mitochondria-derived ROS and apoptosis [66,67]. 
In summary, the observed changes in the abundance of 
glycolytic proteins may be a cold-induced adaptation 
that promotes NADPH reduction while maintaining 
ATP generation and minimizing mitochondria-generated 
oxidative stress. 

Proteins that regulate gene expression 

Cold-adaptation induced significant changes in the 
abundance of proteins that regulate gene expression. 
Histones and histone modifications, such as methylation 
or acetylation, play an important role in regulating tran- 
scription by modifying chromatin structure [68]. The 
histone H2A.V increased significantly at 25°C. In Arabi- 
dopsis and in budding yeast cells, incorporation of the 
closely related variant H2A.Z [69] into nucleosomes is a 
requirement for direct temperature-induced changes in 
gene expression [70]. Whether or not H2A.V plays a 
similar role in promoting the expression of cold-respon- 
sive endothelial genes is unknown. 

Histone methylation mediates short-term chromatin 
condensation and transcriptional repression (DNA 
methylation affects the long-term silencing of genes) 
[71]. In addition to its role in generating cysteine for 
GSH synthesis, S-adenosylhomocysteine hydrolase 
(SAHH) is critical to histone/DNA methylation and 
gene silencing because it eliminates S-adenosylhomocys- 
teine (SAH), a by-product and potent inhibitor of S-ade- 
nosylmethionine (SAM)-dependent methyltransferase 
reactions [72] (SAM is the methyl donor). The increase 
in SAHH observed at 25°C may therefore promote his- 
tone/DNA methylation and potentially gene silencing 
[73]. There is also experimental evidence that increases 
in Nampt, which we observed at 25°C, promotes histone 
deacetylation and consequently gene silencing [74]. The 
level of NAD and Nampt in mouse fibroblasts was 
shown to regulate the activity of Sir2, a protein deacety- 
lase [31]. Similarly, the overexpression of Nampt in 
smooth muscle cells enhanced cell resistance to oxida- 
tive stress and increased the activity of SIRT1, the 
human homolog of Sir2 [75]. SIRT1 regulates nucleo- 
some and chromatin structure by deacetylating histones 
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H3 and H4 and mediates stress resistance, apoptosis and 
inflammatory responses by deacetylating the transcrip- 
tion factors FOXO, p53 and NF-kB [76]. The increase 
in the expression of SAHH and Nampt at 25°C, and 
their capacities to modulate both transcriptional repres- 
sion and resistance to oxidative stress, suggests that 
these proteins may have pivotal roles in the adaptation 
to cold. 

Additional findings point to an overall decrease in tran- 
scription and translation at 25°C. There were significant 
decreases in the abundance of proteins associated with 
transcription initiation, mRNA elongation, processing 



and nonsense- mediated decay, and stress granule assem- 
bly. The inhibition of IMP dehydrogenase has been 
shown to disrupt DNA and RNA synthesis [42] and its 
decrease at 25°C potentially limits global transcription 
and translation. The increase in transcriptional repressors 
ENOl and HNRNPDL at 25°C complement earlier find- 
ings that showed increases in the repressors prohibitin, 
L3MBTL2 and Zinc finger protein 224 [52]. Although 
ENOl (a-enolase) is recognized as a glycolytic protein, 
an alternate translation product of the ENOl mRNA is a 
repressor of the c-myc gene [77]. The myc family of pro- 
teins are transcription factors that regulate ribosome 
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biogenesis, protein synthesis, cell growth and prolifera- 
tion [78]. Although we did not identify c-myc, N-myc 
was identified and it had the largest fold-change decrease 
(-2.38) of any protein detected in this study (Additional 
file 1; N-myc was excluded from the analysis because of 
its classification as a priority 4 protein [24,26]). However, 
the decrease in N-myc abundance and increase in ENOl 
match the apparent decrease in ribonucleoprotein com- 
plex biogenesis (11 out of 11 proteins decreased), transla- 
tion (38 out of 40 decreased) and unfolded protein 
binding (9 out of 9 decreased). By contrast, cold-adapta- 
tion increased the abundance of RBM3, a cold-inducible 
RNA-binding protein that increases the efficiency of glo- 
bal protein synthesis at 37°C and at mild hypothermic 
temperatures [42]. It has been suggested that RBM3 may 
prevent a disproportionate cold-induced decrease in 
translation [42]. The fact that one third of the changes in 
protein abundance at 25°C were increases suggests that 
RBM3 may have an important role in facilitating these 
increases. 

Response to stress proteins 

The response to stress proteins (Figure 6) that increased or 
decreased significantly with cold-adaptation are associated 
with processes like chromatin binding, transcription, DNA 
repair, pre-mRNA processing, regulation of translation, 
protein folding and cell adhesion. Five of the identified 
proteins are found in the ER, which is the site where secre- 
tory and cell surface proteins are synthesized, folded and 
modified [79] . The decrease in the abundance of these ER- 
resident proteins after 72 h at 25°C is additional evidence 
that protein synthesis and protein folding is repressed in 
the cold-adapted state. Excess unfolded or misfolded pro- 
teins in the ER activates the Unfolded Protein Response or 
UPR. The UPR initiates translational repression, increases 
the transcription of ER chaperone genes and increases the 
degradation of misfolded proteins to re-establish the equi- 
librium between protein synthesis and protein folding 
reactions [79]. It is not known if ER function is impaired 
at 25°C or if translational repression due to hypothermia is 
initiated by the UPR, but evidence of UPR activation in 
the tissues of hibernating mammals has been reported 
[80]. Activation of the UPR is also a potential mechanism 
for increasing GSH in cold-adapted cells [16] because of 
the important role that GSH plays in regulating redox 
homeostasis during the oxidative folding of proteins [29] . 
ER stress-activated transcription factors ATF4 and Nrf2 
increase GSH synthesis [81,82] and increases in endothe- 
lial GSH have been linked to hyperoxic- or cycloheximide- 
induced suppression of translation [83]. 

YB-1 and CSDA are two nucleic acid-binding stress 
proteins that contain a central cold-shock domain, a 100- 
amino acid sequence that derives its name from its 
resemblance to bacterial cold-shock proteins [38]. In 



homeothermic organisms, cold-shock domain (CSD) 
proteins have been associated with the transcriptional 
regulation of stress-response, growth factor, and prolif- 
eration-associated genes [38]. In bacteria, cold-shock 
proteins function as mRNA chaperones that facilitate 
translation [84]. Their expression levels increase during 
the acclimation to cold and then later decline below pre- 
cold-shock levels [85]. It is this archaic function of CSD 
proteins that may have special relevance to the endothe- 
lial adaptation to cold. 

Cold-adapted cells had significantly higher levels of 14- 
3-3 proteins zeta/delta (Figure 6, 7) and theta (Figure 5). 
14-3-3's are a family of conserved proteins that control 
many diverse processes by binding to specific phosphory- 
lated sites on target proteins to induce changes in the tar- 
get conformation or to change the target's interactions 
with other proteins [86]. In plants, 14-3-3's help regulate 
processes such as photosynthesis, ATP production and 
nitrate reduction in response to changing environmental 
factors including cold .stress [87]. In mammalian cells, 14- 
3-3's regulate responses to growth factors, including the 
stimulation of glycolysis, inhibition of apoptosis and inhi- 
bition of transcription [86] . 

Cytoskeletal proteins 

The adaptation of endothelial cells to cold induced signifi- 
cant changes in the abundance of cytoskeletal proteins 
that implicate diminished microtubular and actin net- 
works, an expanded intermediate filament network and 
modified cell-cell and cell-extracellular matrix interactions. 
The actin cytoskeleton underlies the plasma membrane to 
provide support, resist tension and maintain cell shape. 
However, it is a highly dynamic structure that responds to 
environmental stimuli and rapidly reorganizes to change 
cell shape, form membrane protrusions, (such as lamelli- 
podia and filopodia for cellular locomotion), or form stress 
fibers (bundles of radially oriented actin-filaments cross- 
linked with myosin motor proteins) for retracting the 
plasma membrane. The proteomic findings suggest that 
there is a significant reorganization of the actin cytoskele- 
ton in response to the reduced temperature. Decreases in 
actin nucleation and cross-linking proteins (ARPC3 [43] 
and FLNA [44]), increases in actin filament severing pro- 
teins (CFL1 and DSTN [45]) and in proteins that bind to 
G-actin to inhibit polymerization (CAP1 [46] and 
TMSB10 [47]), allude to a net disassembly of the actin 
cytoskeleton at 25°C. Disassembly is further indicated by 
the increase in RhoGDI2, which promotes actin depoly- 
merization [49], and by the decrease in septin 2, whose 
depletion leads to dissociation of actin stress-fibers [48] . 

Microtubules consist of a- and |3 -tubulin heterodimers 
and form robust structures that are resistant to bending 
or compression [88]. Microtubules organize the cyto- 
plasm, position the nucleus and organelles, transport 
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organelles and make up the mitotic spindle [88]. The 
decrease in a-tubulins, microtubule associated proteins 
and chaperonins at 25°C indicate a decrease in the micro- 
tubule network (11 out of 12 microtubule-associated pro- 
teins decreased). Our preliminary studies showed a 
matching increase in the microtubule stabilizing protein 
LIS1 at 25°C [52], which points to a decrease in microtu- 
bule catastrophe events and diminished microtubule 
dynamics [89]. One consequence of these changes may 
be a reduction in GTP utilization because microtubule 
assembly/disassembly requires the hydrolysis of tubulin- 
bound GTP to GDP [88]. Such a reduction would match 
the apparent decrease in translation, which is another 
GTP-utilizing process, and the decrease in IMPDH2 (Fig- 
ure 2), which is the rate-limiting enzyme in de novo 
GMP synthesis. 

In contrast to diminished microtubular and actin-fila- 
ment networks, there was a cold-induced increase in type 
III intermediate filaments (IFs; Figure 7) that form a third 
class of cytoskeletal network. Cells at 25°C developed a 
significant increase in vimentin, a major IF that provides 
mechanical support for the cell and takes part in addi- 
tional processes such as the regulation of cell adhesion to 
the extracellular matrix. A second IF, peripherin, also 
increased. IFs are highly dynamic structures that self- 
assemble without expenditure of GTP or ATP required by 
tubulin or actin polymerization, respectively. It has been 
theorized that since actin dynamics may account for up to 
50% of ATP turnover in resting cells, a change in cytoske- 
letal structure to a less dynamic form may reflect a 
decrease in the energy state of the cell or possibly repre- 
sent a form of ATP preservation for more crucial activities 
[90]. The substitution of an expanded IF network for 
diminished actin filament and microtubule networks may 
therefore be a possible sign that cold-adapted cells are in 
an energy-reduced or energy-conserved state. 

Cytoskeletal changes such as those noted above are fre- 
quently mediated through changes in cell-cell and cell- 
matrix adhesion and our study indicates that similar pro- 
cesses are occurring with cold-adaptation. Cdc42, a Rho 
GTPase that regulates actin cytoskeletal remodeling in 
response to environmental signals, increased significantly 
at 25°C. Activation of Cdc42 preserves the integrity of 
endothelial cell-cell junctions [91], prevents membrane 
retraction [91] and induces filopodia formation [92]. 
These properties may contribute to the improved ability of 
cold-adapted endothelial cells to reverse intercellular gap 
formation during rewarming from 0°C [16]. The increase 
in Cdc42 in our study coincided with significant increases 
in fascin, a protein that organizes actin filaments into the 
parallel bundles that are incorporated into filopodial pro- 
trusions [93]; in myosin light chain subunits (MYL6, 
MYL6B, MYL12A and MYL12B) that regulate actinomyo- 
sin contractility and filopodia formation [94]; and in 



fibronectin 1 and thrombospondin 1 (Figure 6), which are 
secreted components of the extracellular matrix. The asso- 
ciation of fascin with actin is regulated by the binding of 
integrin, a transmembrane receptor, to the extracellular 
matrix and to thrombospondin and fibronectin in particu- 
lar [95]. In addition, there were changes in the abundance 
of proteins that link integrin to the cytoskeleton. For 
example, CasL, which mediates signals between integrin 
and possibly vimentin [96], and talin 1, which links integ- 
rin to actin [97], increased significantly. By contrast, plec- 
tin, a linker protein between IFs and microtubules, 
microfilaments or integrins [98], and kindlin-3, which is a 
co-activator of integrin [99], decreased at 25°C. Cold- adap- 
tation also caused changes in components of cell-cell com- 
plexes. There were significant decreases in |3-catenin, 
which links the adherens cell-cell junction to F-actin 
[100]. Our previous findings showed decreases in the 
membrane receptor EphA3, which is a member of the Eph 
receptor family of receptor tyrosine kinases that induce 
changes in cell shape and movement upon the binding of 
ephrin ligands from an adjacent cell [101]. Overall, these 
findings suggest that cytoskeletal remodeling during cold- 
adaptation may be stimulated by increases in fibronectin 
and thrombospondin secretion, integrin binding and 
cdc42 activation, and by decreases in adherens junction 
and ephrin/ephrin receptor cell-cell signaling. 

Lipids/lipid fluidity 

Organisms that successfully adapt to cold change the 
lipid composition of their cell membranes to maintain 
membrane fluidity and function with decreasing tem- 
perature. Some cold-induced membrane changes 
observed in fish are an increase in the proportion of 
unsaturated to saturated fatty acids, an increase in fatty 
acid chain length and an increase in the proportion of 
phosphatidylethanolamine to phosphatidylcholine [102]. 
These changes are brought about by adjustments in lipid 
metabolism, including de novo synthesis, elongation and 
desaturation of fatty acids; de novo phospholipid synth- 
esis and desaturation; and phospholipid catabolism and 
turnover in the membrane [102]. In our study, few lipid 
metabolic proteins were identified so there is little infor- 
mation about potential modifications to endothelial cell 
membranes in response to cooling. Fatty acid synthase 
was significantly down regulated in cold-adapted cells 
(Figure 2), which suggests that de novo fatty acid synth- 
esis is not essential. Interestingly, there was a cold- 
induced increase in phospholipase A 2 receptor 1 
(PLA2R1 was a priority 4 protein; see Additional file 1), 
which binds the lipolytic enzyme phospholipase A 2 
(PLA 2 ). PLA 2 s hydrolyze phospholipids at the sn-2 posi- 
tion, releasing free fatty acids and lysophospholipids from 
the membrane. The sn-2 position is frequently occupied 
by polyunsaturated fatty acids (PUFA), so PLA 2 activity 
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potentially regulates membrane unsaturation and fluidity. 
In support of this hypothesis, Zhang et al. demonstrated 
that inhibition of Ca 2+ -independent PLA 2 activity 
increased PUFA in membrane phosphatidylcholines [103]. 
PLA 2 activity at hypothermic temperatures is well-known 
[104-106]. It is not known, however, if PLA2R1 has a sti- 
mulatory or inhibitory effect on endothelial PLA 2 s. 

In addition to a direct effect of cooling on membrane 
fluidity, there may be an indirect effect that is caused by 
cold-induced oxidative stress and peroxidation of unsa- 
turated lipids. The lipid peroxidation products 4-hydro- 
xynonenal and malondialdehyde have been shown to 
increase the rigidity of cell membranes [107,108] and 
have been measured following cold exposure and 
rewarming [50,51]. The increase in peroxiredoxin, alde- 
hyde dehydrogenase and glutathione S-transferase in 
cold-adapted cells potentially attenuate lipid peroxida- 
tion in endothelial membranes and may help preserve 
membrane fluidity during cooling. 

Comparisons to earlier proteomic studies 

This study was performed using a well-validated label- 
free mass spectrometry-based method [24,26,109,110] 
that is rapid, sensitive and has a three- to four-fold 
higher protein dynamic range (protein mass and abun- 
dance) than two-dimensional gel electrophoresis (2-DE) 
[111]. Almost one third of the 1089 proteins that were 
identified significantly increased or decreased in abun- 
dance during cold-adaptation and 181 proteins of high 
identification confidence were included in the functional 
analysis. This number greatly exceeds the number of 
proteins identified and analyzed in previous 2-DE-based 
hypothermia studies [112,113]. Ten of 26 proteins iden- 
tified by Baik et al. [112] changed significantly in our 
study, including the glycolytic proteins pyruvate kinase 
M2, phosphoglycerate kinase and glyceraldehyde-3- 
phosphate dehydrogenase; cytoskeletal proteins a-tubu- 
lin, annexin Al and vimentin; chaperones HSP90|3 and 
HSPA8; and the ER-resident protein disulfide isomerase. 
With the exception of PDI and HSPA8, the changes in 
abundance were in the same direction as in our study. 
Kumar et al. [113] observed similar cold-induced 
increases in the abundance of vimentin, aldehyde dehy- 
drogenase 1A1 and GAPDH as were reported in this 
study. They also reported decreases in tubulins (P -tubu- 
lin rather than a-tubulin) and annexins (A4 rather than 
Al) [113]. In addition, Kumar et al. [113] identified a 
decrease in a-enolase (we noted an increase) and an 
increase in the cytoskeletal chaperone CCT2 (we 
reported decreases in CCT4 and 5). 

Conclusions 

Differential proteomic analysis of endothelial cells col- 
lected before or after cold-adaptation at 25°C 



demonstrated significant changes to a broad range of 
biological processes and a new phenotype of enhanced 
redox protection and down-regulated anabolic and 
energy-consuming processes. The principal findings are 
depicted in Figure 9 and provide a preliminary model of 
the endothelial adaptive response to cold and a frame- 
work to guide future investigations. Cold-adaptation 
increases the abundance of redox active proteins thiore- 
doxin, thioredoxin reductase, glutathione S-transferase 
and tyrosine hydroxylase [52] and decreases protein dis- 
ulfide isomerases 1 and 4 and glutaredoxin-3 [52]. 
These changes correlate with a recovery of protein 
thiols, an increase in protein glutathionylation and 
diminished oxidative stress when cells are rewarmed 
from 0°C. Protein glutathionylation may prevent irrever- 
sible oxidation and formation of intermolecular disulfide 
bonds and therefore represents a potentially important 
protective function of GSH in cold-adapted cells. The 
increase in abundance of S-adenosylhomocysteine 
hydrolase (SAHH), a component of the methionine- 
cysteine transsulfuration pathway that generates cysteine 
for GSH synthesis, may play a role in generating the 
higher concentrations of GSH observed in cold-adapted 
cells [16]. In addition to protecting proteins from oxida- 
tive injury, cold-adapted cells may be more resistant to 
lipid peroxidation than non-adapted cells because of 
increases in peroxiredoxin 1 and aldehyde dehydrogen- 
ase (ALDH), which enzymatically reduce lipid peroxides 
and detoxify aldehydes, respectively. The increase in 
antioxidant activities requires the availability of reducing 
equivalents from NADPH and the proteomic data show 
several adaptive changes that could increase NADPH 
availability. Cold-adapted cells have higher levels of 
Nampt, the rate-limiting enzyme in NAD synthesis from 
nicotinamide. A significant decrease in PFK, the rate- 
limiting enzyme of glycolysis, potentially redirects glu- 
cose catabolism to the pentose phosphate pathway to 
increase NADP + reduction to NADPH. The increase in 
ALDH abundance (and presumably activity) also poten- 
tially increases NADP + reduction. Most importantly, the 
attenuation of NADPH-dependent anabolic processes 
such as DNA, RNA and fatty acid synthesis due to 
decreases in RRM1, IMPDH2, and FAS expression, 
respectively, could further increase NADPH availability 
for antioxidant activities. 

Central to the proteomic changes is the evidence of a 
regulated decrease in energy consuming processes such 
as protein synthesis, DNA synthesis, fatty acid synthesis 
and cytoskeletal dynamics. The net decrease in proteins 
associated with translation, ribonucleoprotein complex 
biogenesis, RNA processing, mRNA metabolism and 
protein folding at 25°C is consistent with a generalized 
cold-induced inhibition of protein synthesis described 
by others [19]. Increases in transcriptional repressors 
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Figure 9 Endothelial adaptive response to moderate hypothermia (25°C/72 h). See Conclusions for details. Thioredoxin (Trx), thioredoxin 
reductase (TrxR), tyrosine hydroxylase (TH), reduced glutathione (GSH), glutathione S-transferase (GST), protein disulfide isomerase (PDIA1 and 4), 
reduced protein thiols (P(SH) 2 ), protein disulfide (PSSP), glutathionylated protein (P-SSG), oxidized protein (P-S-), S-adenosylhomocysteine 
hydrolase (SAHH), peroxiredoxin 1 (Prx1), aldehyde dehydrogenase (ALDH), nicotinamide phosphoribosyltransferase (Nampt), reduced/oxidized 
nicotinamide adenine dinucleotide phosphate (NADPH/NADP + ), fatty acid synthase (FAS), ribonucleotide reductase M1 (RRM1), inosine-5'- 
monophosphate dehydrogenase 2 (IMPDH2), phosphofructokinase (PFK), pentose phosphate pathway (PPP), histone (H2AV, H2BFC, H2BFQ), 
histone binding proteins (NASP or nuclear autoantigenic sperm protein), chromatin binding proteins (High-mobility group protein B1 (HMGB1), 
high-mobility group protein 1 - 1 i ke 1 0 (HMG1 L10), barrier to autointegration factor 1 (BANF1)), RNA binding motif protein 3 (RBM3), Y-box binding 
protein 1 (YB-1), cold-shock domain protein A (CSDA), prohibitin (PHB), enolase 1 (EN01), lethal(3)malignant brain tumor-like 2 protein (L3MBTL2), 
heterogeneous nuclear ribonucleoprotein D-like (HNRNPDL), Zinc finger protein 224 (ZNF224). 



and changes in the abundance of histone and non-his- 
tone chromatin-binding proteins are possible indicators 
of a concurrent decrease in transcription. Changes in 
the expression of cold-shock proteins (RBM3 increased, 
CSDA and YB-1 decreased) and the increase in histone 
H2A.V at 25°C potentially counteract the generalized 
decrease in gene expression by enhancing translation 
efficiency and the direct responsiveness of genes to cold, 
respectively. This may be an important requirement for 
increasing the synthesis of a subset of proteins at 
reduced temperatures- one third of identified proteins 
increased in abundance in our study. Two enzymes we 
have identified, SAHH and Nampt, may play pivotal 
roles in mediating adaptation to cold because each 
induces transcriptional repression and increases the 



cell's reducing capacity. SAHH facilitates SAM-depen- 
dent DNA and histone methylation and consequently 
gene silencing [73] while up-regulating GSH synthesis. 
Nampt regulates histone deacetylation and transcrip- 
tional repression [74] while increasing NAD synthesis. 
The ability to down-regulate cellular energy demand 
while increasing the cell's reducing capacity would be 
clinically beneficial where the application of cold-ische- 
mia would otherwise create energy insufficiency, uncou- 
ple metabolic integration and generate oxidative stress 
and injury. Whether or not a similar adaptive response 
occurs in other cell types, such as myocardial or neural 
cells, or following an ischemic insult, such as during the 
application of therapeutic hypothermia, remain impor- 
tant questions. 
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Methods 

Cell culture 

Proliferating human coronary artery endothelial cells (Clo- 
netics; passages 4 to 8) were seeded at a density of 1.4 x 
10 4 viable cells/cm 2 (6-well culture plates, 12.5 cm 2 , 
25 cm 2 or 75 cm 2 culture flasks; Falcon, Franklin Lakes, 
New Jersey) and were maintained at 37°C (95% air + 5% 
CO2) in culture medium (EGM2-MV; Clonetics). Cells 
formed a confluent monolayer after 3 to 4 days of culture 
and were then used in experiments. Cells were adapted to 
cold at 25°C for up to 72 h in a water-jacketed incubator 
(3.8% C0 2 + 96.2% air) and the culture medium was 
replaced with fresh medium every 24 h. Controls consisted 
of confluent cells collected at the start of the adaptation 
period or following continued maintenance at 37°C for up 
to 72 h. For cold storage, cells were washed with Hank's 
balanced salt solution, flasks received endothelial basal 
medium (EBM-2, Clonetics) and were immersed in a cir- 
culating water bath at 0°C. 

Label-free protein quantification method 

A large scale quantitative proteomics experiment was per- 
formed using a well-validated label-free mass spectrome- 
try-based method [24,26,109,110]. Cold-adapted and 
control cells were prepared as previously described [110]. 
Proteins were extracted from the tissue culture cells using 
8 M urea, reduced by triethylphosphine, alkylated by 
iodoethanol, and digested by trypsin [114]. Tryptic pep- 
tides (-20 ug/injection) were analyzed using Thermo- 
Fisher linear ion-trap mass spectrometer (LTQ) coupled 
with a Surveyor HPLC system. C-18 reverse phase column 
(i.d. = 2.1 mm, length = 50 mm) was used to separate pep- 
tides with a flow rate of 200 uL/min. Peptides were eluted 
with a gradient from 5 to 45% acetonitrile developed over 
120 min and data were collected in the triple-play mode 
(MS scan, zoom scan, and MS/MS scan). The acquired 
data were filtered and analyzed by a proprietary algorithm 
that was developed and described by Higgs, et al. [26,109]. 
Database searches against the IPI (International Protein 
Index) human database was carried out using both the X! 
Tandem and SEQUEST algorithms. Protein quantification 
was also carried out using the same proprietary algorithm 
described previously [24,26,110]. 

Bioinformatics analysis 

IPIs of Priority 1 and 2 proteins were converted to Uni- 
Prot Accession numbers using the EMBL-EBI PICR Ser- 
vice, http://www.ebi.ac.uk/Tools/picr/. Accession 
numbers were then uploaded to the DAVID (Database 
for Annotation, Visualization and Integrated Discovery) 
Bioinformatics Resources 6.7 website http://david.abcc. 
ncifcrf.gov/tools.jsp and the proteins were classified and 
organized according to their GO Biological Process, Mole- 



cular Function, Cell Component and SP-PIR Keyword 
annotations using the Gene Functional Annotation Tool. 

Measurement of protein thiols (PSH) and mixed disulfides 
(PSSG) 

Confluent HCAEC from T-75 flasks were rinsed twice 
with HBSS and scraped in 1.5 ml of 5% TCA in 10 mM 
HC1 with 5 mM EDTA. Protein pellets and non-protein 
supernatants were separated by centrifugation and 
stored at -80°C until analysis of protein thiols (PSH) and 
mixed disulfides (PSSG) using the method described by 
Nagy [115] with some modifications. The protein pellets 
were extracted with diethyl ether to eliminate TCA, 
redissolved in 0.25 ml 6 M GuHCl, and pH adjusted to 
8 by addition of 0.5 ml phosphate buffer with 5 mM 
EDTA, pH 8.6. Insoluble material was removed by cen- 
trifugation and supernatants were used to determine the 
levels of PSH and PSSG. For PSH, 100 ul of the sample 
was reacted with 0.9 ml of 0.1 mM 5,5'-dithio-bis-(2- 
nitrobenzoic acid) (DTNB, 1 mM stock in methanol) in 
125 mM sodium phosphate buffer containing 6.25 mM 
EDTA, pH 8.6. Absorbance was measured at 412 nm 
after 30 minutes at room temperature. GSH standards 
(50-500 uM) were similarly measured and PSH-levels in 
samples were calculated directly from the standard 
curve. 

PSSG levels were determined in the reaction mix com- 
prised of 0.1 ml sample, 0.9 ml of 125 mM sodium 
phosphate buffer, pH 7.5 (with 6.25 mM EDTA) con- 
taining 0.1 mM DTNB, 0.03 mM NADPH and 20 ul of 
glutathione reductase (25 units per ml). The absorbance 
was measured at 412 nm after 10 minutes of incubation 
at room temperature. PSSG levels in samples were cal- 
culated from a GSH standard curve (0.5 - 10 uM) run 
under same conditions. 

Glutathione determination 

Confluent cells were washed two times in HBSS and 
then lysed in 500 uL of 125 mM sodium phosphate buf- 
fer with 6.3 mM EDTA and 2% Triton X-100. Protein 
was removed from lysates by precipitation with 500 uL 
of metaphosphoric acid (10% w/v) and centrifugation at 
14,000 x g for 5 minutes. The supernatant was neutra- 
lized with 5 N sodium hydroxide and stored at -80°C 
until t-GSH and GSSG was measured using the method 
of Ishii [116]. To determine t-GSH, 100 uL of the neu- 
tralized supernatant was incubated at 25°C with 800 uL 
of 0.3 mM NADPH in 125 mM sodium phosphate buf- 
fer with 6.3 mM EDTA, 100 uL of 6 mM 5,5'-dithiobis 
(2-nitrobenzoic acid) (DTNB) and 20 uL of 25 U/mL 
GSH reductase. After precisely 15 minutes, the optical 
density was measured at 412 nm. To determine GSSG, 
GSH in samples was derivatized by adding 2.5 uL of 
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4-vinylpyridine (95%) per 250 uL solution, mixing vigor- 
ously for 1 min and incubating samples at room tem- 
perature for 30 min. Samples were then centrifuged at 
14,000 x g for 5 min and GSSG was measured in super- 
natants by the procedure described for t-GSH. The mea- 
sured optical densities were converted to concentrations 
using standard curves prepared with GSH (0.5-10 uM) 
or GSSG (0.25-5 uM). 

Measurements of thioredoxin, thioredoxin reductase, and 
glutathione S-transferase activity 

Confluent HCAECs from T-75 flasks were scraped and 
resuspended in 0.15-0.5 ml of assay buffer and cell 
lysates were prepared by sonicating samples two times 
on ice (10 seconds, amplitude of 40) using an Ultrasonic 
Processor. Insoluble materials were removed by centrifu- 
gation at 13,000 rpm for 20 minutes at 4°C and super- 
natants were used to measure enzyme activities and 
total protein concentrations (Micro Bicinchoninic Assay 
Kit; Pierce Thermoscientific). Thioredoxin (Trx) activity 
in samples was determined using an assay kit (Redoxica, 
Littlerock, AK) and following the manufacturer's recom- 
mendations. Briefly, Trx activity was determined by 
measuring the rate of depletion of NADPH at 340 nm 
in the presence of oxidized insulin and thioredoxin 
reductase (TrxR). Reaction blanks with each sample 
were run without the insulin substrate to determine 
Trx-independent NADPH depletion. Positive controls 
consisting of thioredoxin supplied with the kit were run 
to validate the assay. Reaction rates were calculated 
from the slopes of logarithmic plots of A 340 nm versus 
time (in min). The difference in reaction rates between 
the samples and corresponding reaction blanks repre- 
sented Trx activity and was expressed as nmoles of Trx 
activity/ min/ mg protein. 

TrxR activity was measured using an assay kit (Cay- 
man Chemical Company, Ann Arbor, Michigan, USA). 
The assay is based on measuring the rate of reduction 
of DTNB (5,5'-dithiobis(2 dinitrobenzoic acid)) to a 
colored product by TrxR, in the presence of NADPH. 
TrxR-independent reduction of DTNB was determined 
by adding aurothiomalate (ATM), a specific inhibitor of 
TrxR, to the reaction mixture. TrxR activity was deter- 
mined by measuring the difference in the rate of 
increase in absorbance at 412 nm between samples with 
or without ATM and was expressed in nmoles/min/mg 
proteins, using an extinction coefficient of 6.35 mM 1 . 

GST activities were assessed using an assay kit from 
Cayman and following the manufacturer's protocol. GST 
activity was quantified by measuring the rate of conjuga- 
tion of reduced glutathione by GST to CDNB (1-chloro- 
2,4 dinitrobenzene) at 340 nm and was expressed in 
nmoles of conjugated CDNB/min/mg protein, using an 
extinction coefficient of 0.0053 uM~ . 



Additional material 



Additional file 1: Summary of identified proteins Each protein is listed 
with its Rank, significant change between control and cold-adapted cells 
(Signif; Yes or No), Priority, ProteinJD, Protein Name, fold-change (Fold; 
positive represents an increase in cold-adapted cells; negative represents a 
decrease), % ID confidence of the best sequence (961 D), qValue (comparing 
control to hypothermia-treated), mean of control (mean_C; n = 5), mean of 
cold-adapted (mean_H; n = 5), % Coefficient of Variation (%CV; standard 
deviation/mean), number of peptides quantified (Peptide*) and peptide 
sequence with highest ID confidence (Best_Sequence). A detailed 
description of the parameters are in [24]. 

Additional file 2: Summary of biological process, molecular function 
and cell component categories that changed significantly with cold- 
adaptation. Each GOTERM or SP_PIR_KEYWORD is listed with the protein 
count, percentage of the total (178 proteins), P-value, gene identifiers, 
number of genes in the category, total number of genes and fold- 
enrichment. 
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